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ABSTRACT 


This  report  investigates  the  effectiveness  of  the  minimum- 
power  array  p  ^cessing  technique  in  determining  seismometer  inequaliza¬ 
tions.  The  technique  involves  partitioning  the  seismometer  array  into  two 
groups  and  designing  MCF's  for  each  group  so  that  the  tman-square-error 
between  the  two  MCF  outputs  is  a  minimum  under  the  constraint  that  the  out¬ 
put  power  of  one  of  the  MCF's  is  unity.  The  two  MCF  sets  a^e  known  as  the 
group-coherence  filters;  the  difference  between  these  sets  is  knov/n  as  the 
minimum-power  array  processor. 

Estimates  of  the  noise  wavenumber  spectrum  from  the  wave- 
number  responses  of  the  group-coherence  filters  are  distorted  due  to  seis¬ 
mometer  inequalization;  however,  a  more  reasonable  estimate  of  the  noise 
wavenumber  spectrum  from  the  wavenumber  response  of  the  minimum-power 
array  processor  should  be  possible  because  the  minimas  in  the  processor's 
wavenumber  response  correspond  to  the  wavenumber  regions  where  the  wave- 
number  responses  of  the  two  group-coherence  MCF's  are  very  similar  (e.  g.  , 
at  the  peaks  of  the  noise  wavenumber  power  spectrum). 

Seismometer  inequalization  was  to  be  determined  from  the  ad¬ 
justment  in  weight  and  phase  required  for  each  filter  so  that  the  wavenumber 
responses  of  the  group-coherence  MCF's  would  agree  with  a  reasonable  noise 
wavenumber  spectrum.  However,  results  from  the  TFO  long-noise  sample 
and  two  synthetic  models  show  that  the  technique,  although  excellent  for  gen¬ 
erating  maximum  coherent  channels,  lacks  the  wavenumber  resolution  de¬ 
sired  for  studying  seismometer  inequalization.  This  latter  conclusion  is  at 
least  true  for  small  arrays  such  as  TFO. 
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ACRONYMS 


MCF 

TFO 


Multichannel  Filter 

Tonto  Forest  Seismological  Observatory 
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SECTION  I 
INTRODUCTION 


In  applying  linear  least-mean-square-error  theory  to  a  single¬ 
channel  prediction  problem,  a  reference  channel  is  estimated  by  linearly 
combining  a  group  of  channels  so  that  the  mean- square-error  is  a  minimum. 
The  minimum-pewer  array  theory  is  a  generalization  of  this  concept  where 
the  reference  channel  itself  is  a  linear  combination  of  another  group  of  chan¬ 
nels  x  and  has  unity  power. 

In  this  report,  vectors  x  and  ^  represent  two  ordered  sets  of 
complex  Fourier  transforms  obtained  at  a  given  frequency  from  the  two  groups 
of  seismometer  outputs.  There  are  now  two  sets  of  linear  multichannel  fil¬ 
ters  (MCF's):  one  operates  on  group  x  to  generate  the  reference  channel; 
the  other  operates  on  groups  to  predict  the  reference  channel. 

Minimization  of  the  mean- square-error  by  varying  the  two 
sets  of  MCF's  leads  to  the  solution  of  a  generali^d  eigenvalue  matrix  equa¬ 
tion.  The  minimizing  pair  of  multichannel  filters  are  known  as  the  group- 
coherence  filters.  The  multichannel  filter  formed  from  the  difference  be¬ 
tween  the  two  group- coherence  filters  is  called  the  minimum-power  array 
processor. 

Mean-square-error  is  numerically  equal  to  the  fraction  of  the 
power  in  the  normalized  reference  channel  (the  channel  obtained  by  applying 
the  group-coherence  filter  to  the  set  of  channels  x)  which  cannot  be  linearly 
predicted  from  group  jr.  The  numerical  value  of  the  group  coherence  is  de¬ 
fined  to  be  the  predictable  fraction  of  the  power  in  the  reference  channel; 
thus,  a  minimum  mean-square-error  is  associated  with  a  maximum  group 
coherence,  and  vice  versa. 
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An  important  property  of  group  coherence  is  its  invariance 
with  respect  to  any  nonsingular  linear  transformation  of  the  channels  within 
the  two  groups;  i.e.  ,  group  coherence  is  unchanged  by  scaling,  by  frequency 
filtering,  or  by  combining  the  channels  within  a  group  by  any  linear  rever¬ 
sible  network  filter.  For  seismometer  arrays,  the  group  coherence  between 
two  arrays  is  independent  of  any  inequalization  problems  possessed  by  the 
seismometers  in  the  arrays;  however,  if  seismometer  inequalization  is 
severe,  the  wavenumber  responses  of  the  group-coherence  filters  differ  from 
reasonable  wavenumber  power  responses.  Except  for  this  inequalization  ef¬ 
fect,  the  wavenumber  responses  of  the  two  group-coherence  filters  should 
tend  to  peak  and  be  highly  similar  in  regions  where  the  wavenumber  power 
spectrum  of  the  array  data  is  a  maximum.  The  wavenumber  response  of  the 
minimum-power  array  processor,  which  is  the  difference  between  the  two 
group-coherence  filters,  should  have  a  small  power  response  at  the  wavenum¬ 
ber  peaks.  Thus,  highly  coherent  energy  such  as  that  generated  by  storms  or 
earthquakes  would  appear  as  deep  troughs  in  the  wavenumber  response  of  the 
minimum- power  array  processor.  A  reasonable  estimate  of  the  wavenumber 
power  spectrum  should  be  possible  from  the  wavenumber  response  of  the 
minimum-power  array  processor. 

If  the  weights  of  the  group-coherence  filters  were  adjusted  to 
compensate  for  seismometer  inequalizations,  the  filters'  wavenumber  re¬ 
sponses  would  agree  with  the  noise  wavenumber  power  spectrum.  Adjurt- 
ments  could  be  made  by  fixing  the  filter  weight  for  one  seismometer  (the 
center  one,  for  example)  as  unity  and  varying  the  filter  weights  correspond¬ 
ing  to  other  seismometers  so  that  the  difference  between  the  group- coherence 
filter  wavenumber  responses  and  the  estimated  noise  wavenumber  spectrum 
would  be  a  minimum.  The  possibility  that  this  processing  scheme  would  lead 
to  determining  the  amount  of  seismometer  inequalization  motivated  the  in¬ 
vestigation  covered  in  this  report. 
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The  broader  objectives  of  this  study  are 

•  To  investigate  whether  the  wavenumber 
response  of  the  minimum-power  array 
processor  can  be  used  for  estimating  the 
noise  wavenumber  power  spectrum  by  de¬ 
tecting  and  isolating  regions  of  highly  co¬ 
herent  energy 

•  To  investigate  whether  the  group- coherence 
filters  and  their  wavenumber  power  responses 
can  be  used  for  determining  the  amplitude  and 
phase- response  inequalizations  of  the  seis¬ 
mometers  in  the  array 

Data  used  in  the  study  are  from  the  TFO  long-noise  sample 
which  was  the  subject  of  Array  Research  Special  Report  No.  23.  ^  Al6o  used 
are  synthetic  data  modeled  to  resemble  the  TFO  data  but  having  no  seismom¬ 
eter  equalization  problem.  The  following  results  concerning  group  coherence 
have  been  derived  from  this  research. 

The  group-coherence  technique  is  excellent  for  measuring  the 
basic  similarity  between  the  two  arrays  of  seismometers,  and  the  group- 
coherence  filters  are  useful  if  the  objective  is  to  generate  the  maximum 
coherent  MCF  outputs.  The  group- coherence  concepi.  is  of  little  value  in 
determining  seismometer  inequalization  —  at  least,  between  array  groups 
having  small  array  separation,  such  as  TFO.  This  is  because  the  wave- 
number  response  of  the  minimum-power  filter  lacks  the  resolution  needed 
for  estimating  the  noise  wavenumber  power  spectrum. 
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SECTION  II 

MINIMUM- POWER  ARRAY  PROCESSING  THEORY 

A.  GENERALIZED  MULTICHANNEL  PREDICTION 

The  term  "channel"  in  this  report  denotes  a  complex  random 
variable.  The  Fourier  transform  of  the  output  of  a  -seismometer,  since  it 
is  a  complex  random  variable,  is  called  a  seismic  channel. 

Let  an  array  of  m  +  n  seismic  channels  (Figure  II- 1)  be  par¬ 
titioned  into  two  ordered  groups  x  =  jx^  |i  =  1,  ...»  and  =  jyj  |j  =  1»  •  •  •  n| 
where  x  and  x  form  two  complex  multivariate  random  column  vectors.  In 
the  generalized  multichannel  prediction  problem,  group  x  is  linearly  com¬ 
bined  to  form  a  reference  channel  of  unity  power,  and  a  linear  combination 
of  groups  is  used  for  obtaining  a  least-mean-square-error  estimate  of  the 
normalized  reference  channel. 

Let  the  two  arbitrary  sets  of  complex  linear  multichannel  fil¬ 
ters  which  operate  on  groups  x  and  x  be  &  s  fei  U  =  •  •  •  »  mf  and  h  =  ihj  I 

j  =  1,  .  .  .  ,  n},  respectively,  and  let  e  be  the  difference  between  the  outputs 
of  £  and  h.  Note  that  MCF  £  generates  the  reference  channel,  that  MCF  h  is 
the  prediction  filter,  and  that  G  is  a  complex  random  variable  which,  in  the 
future,  will  be  referred  to  as  error.  Thus, 


T  .T 
£  x  -  h  x  =  e 


where  T  is  used  to  denote  conjugate  transpose. 


(2-1) 
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B.  MINIMIZATION  OF  MEAN-SQUARE- ERROR 
The  mean-square-error  is  given  by 


where  a  line  indicates  the  average  value.  Therefore, 


where 


is  autopower  matrix  for  group  x 

is  autopower  matrix  for  grouper 

is  crosspower  matrix  for  groups  x  and 
T 

n  =  q 

“yx  — xy 


2 

It  is  desired  to  find  £  and  h,  which  will  minimize  e  |  ,  under 


the  constraint  that  g  xl  is  unity. 


Keeping  £  fixed,  the  variation  of  U  |  with  respect  to  h  in 
Equation  2-2  gives 


U-3 
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where  0  is  the  null  vector  and  6  denotes  the  variation.  Note  that 


n  n 

=XX  amxy 

n  n 

»yx  =yy 


n  n 

■ocx  *sxy 

n  n 

-yx  -yy_ 


=  -6h  (n  g  -  fl  h^ 

-  V-yx  *  -yy  -/ 


(2-4) 


From  Equations  2-3  and  2-4,  it  is  seen  that 


6 1®  1 2  =  -2  [«iT  (a  £  -  a  *)]  <2-s) 


The  condition  for  |e  |  to  be  stationary  is  that,  for  an  infini- 

T  2  T 

tesimal  change  6h  ,  6  |  e  |  =0.  This  implies  that  the  coefficient  of  6h  in 

Equation  2-5  must  vanish;  i.  e.  , 


n  g  -  n  h  =  0 

— yx  ■“  =yy  — 


n  h  =  n  g 

=yy  —  =yx  ■a 


(2-6) 


Since  n  is  a  covariance  matrix,  it  is  positive  definite  and, 

~yy  . 

hence,  nonsingular;  i.e.  ,  fi^  1  exists.  Premultiplying  both  sides  of  Equa¬ 
tion  2-6  gives 


h  =  n  n  g 

—  =yy  =yx  a 


(2-7) 


H-4 
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Substituting  h  from  Equation  2-7  into  Equation  2-2  and  factor¬ 


ing  gives 


n  -1  ci 

=yy  =yx 


n  ci 

=XX  srxy 


ci  n 

=yx  =yy 


-1 

■n  n 

=yy  =yx 


Sl\T  (ci  -  n  ci  _1  ci  \ 

\  l\=™  =xy  =yy  =yx/ 


■n  q 

=yy  ==yx 


n  n  n  \  g  = 

=xy  =yy  =yx  J  & 


(2-8) 


where  I_  is  the  m  x  m  identity  matrix  and  Q  is  the  m  x  n  null  matrix. 

From  Equation  2-8,  it  is  seen  that  mean-square-error  J g  |  ^ 
depends  only  on  the  choice  of  the  MCF  vector  £  which  generates  the  reference 
channel.  The  problem  now  is  to  find  a  minimizing  MCF  vector  g  which  is 
normalized  to  unity  output  power.  Since  the  matrix  Cl  is  positive  definite 
and  since  ,  the  matrix  «yy"  1  £yX  is  Positive  definite.  The 

matrix=xx’  because  it:  is  a  covariance  matrix,  also  is  positive  definite. 

From  matrix  theory,  if  either  of  the  two  matrices  0  or  Q  q  n 

=xx  =xy  =yy  =yx 

is  positive  definite,  there  exists  a  nonsingular  matrix  which  will  simultanerusly 
reduce  both  the  and  1  Q.yx  matrices  to  a  diagonal  form.  Let 

G  be  such  an  m  x  m  matrix  normalized  so  that 


G  n  G  =  I 
—  =sx  =  = 


(2-9) 
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and 


gt  fn  n  _1  o  1  g  =  |x|' 

*  L-xy  — yy  =yx J  —  1 


(2-10) 


f  r 

where  \  is  a  diagonal  m  x  m  matrix. 


Equations  2-9  and  2-10  are  generalized  eigenvalue  matrix 


equations. 


Since  G  is  an  m  x  m  nonsingular  matrix,  its  column  vectors 
form  a  basis  for  generating  arbitrary  MCF's  for  group  x.  Let  £  be  such  an 
MCF  given  by 


£  =  G  £ 

where  £  =  jc.  |i  =  1,  .  .  .  ,m|  is  a  set  of  arbitrary  complex  scalars.  Normal¬ 
izing  £  to  unity  output  power  gives 


£.  Q  &  =  1 

=xx  “ 


(2-11) 


T  T 

c  G  n  G  c  =  1 
—  =  =xx  =  — 


(2-12) 


simplifies  to 


Since,  from  Equation  2-9,  G  ft  G  =  _I ,  Equation  2-12 

“  “XX  “  — 


T  , 
c  c  =  1 


m  2 

E  S  =  1 


(2-13) 


Equation  2-13  gives  the  condition  for  unity  normalization. 
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- 1  T 

Pre-  and  postmultiplying  £  Qyy  Q  by  £  and  £,  re¬ 


spectively,  yields 


T  - 1  T  T  -1 

&QQ  ng  =  cGnn  n  gc  (2-14) 

=xy  =yy  =yx  —  =  =xy  =yy  =yx  =  — 

T  - 1  2 

Since,  from  Equation  2-10,  G  Q  Q  Q  G  =  Ixl  , 

=  =xy  =yy  — yx  =  1  =  1 


Equation  2-14  reduces  to 


T  „ 

&  Q. 


m 


n 

=xy  =yy 


1fi„l-£T|il2£*  £ 


HI' 


(2-15) 


Subtracting  Equation  2-15  from  Equation  2-11  and  factoring, 

*T  (a** -  axy  C,;1  Qyx)a  - 1  -  f  |c,i2  |x.|2  (2-16, 

Comparing  Equations  2-8  and  2-16,  it  can  be  seen  that 

- 7  m  ?  9 

ur  =  1  -  L  |c.r  i\r  (2-17) 


Equation  2-17  gives  the  mean-square-error  associated  with  an  arbitrary 
normalized  filter  £  =  G  £.  |e  |2  is  minimized  by  setting  c.  =  1  when  |  X.  |2  = 
l\naxl  and  setting  a11  other  c.'s  =  0  in  Equation  2-17;  i.  e.  , 


mm 


=  1  -  1 X  | 

'  max  1 


(2-18) 


C.  GROUP  COHERENCE 


Section  I  states  that  group  coherence  is  measured  by  the  maxi¬ 
mum  predictable  fraction  of  power  in  the  normalized  reference  channel.  Since 

the  reference  channel  is  normalized  to  unity  power  and  since  |e  .  I  is  the 

1  min  1 

minimum  prediction  error,  Equation  2-18  shows  that  the  group  coherence 
equals  the  largest  eigenvalue  I  ^max  I* 
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The  eigenvector  associated  with  the  largest  eigenvalue  is  ob¬ 
tained  by  solving  Equations  2-9  and  2-10,  using  an  iterative  scheme  called 
the  power  method.  (This  scheme  has  been  described  in  Large-Array  Signal 
and  Noise  Analysis  Special  Scientific  Report  No.  13.  If  £m„x  is  such  an 
eigenvector,  then,  from  Equations  2-9  and  2-10, 


^max  =xx  ^max  * 


(2-19) 


fcl  [~n  n  " 1  q  1  gi 

“max  wixy  — yy  —yx  “m 


=  X 

ax  1  max 


(2-20) 


The  eigenvector  JSmax  is  known  as  the  group-coherence  MCF  associated  with 


group  x. 


The  group-coherence  MCF  h  ^  associated  with  group  £  is 


obtained  by  substituting  g^^  into  Equation  2-7: 


h  =  n  n  g 

—max  =yy  =yx  “max 


(2-21) 


The  minimum-power  array  processor  output  is  obtained  by  replacing  £  and 

h  in  Equation  2-1  by  g _  and  h  ,  respectively. 

—  “max  —max 
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D.  AUTO-  AND  CROSSPOWER  SPECTRA  OF  THE  TWO  MAXIMUM- 
COHERENT  CHANNELS 

The  autopower  spectrum  of  the  reference  channel  is  given  by 

Equation  2-19.  To  obtain  the  autopower  spectrum  of  the  estimated  reference 

channel,  take  the  conjugate  transpose  of  Equation  2-21  while  recalling  that 

fl  -1  .  „  .  .  T 

-yy  is  Hermitian  and  fi  =Q  .  Thus, 

=yx  =xy 

—max  ”  ^ax  Qxy  =yy  (2-22) 

By  postmulti plying  Equation  2-18  on  both  sides  by  Q  h 
,  ,  -1  =yy  —max 

and  noting  that  £yy  =  I ,  where  I'isannxn  identity  matrix, 

T  T 

“max  =yy  ^max  "  ^max  =xy  \tiax  (2“23) 

By  replacing  h^^  on  the  right-hand  side  of  Equation  2-23  by 
its  equivalent  expression  from  Equation  2-21, 

X  T  •  1 

^max  Oyy  —max  "  &max  £xy  £yy  Syx  4nax 

^maxT  (=xy  =yy  =*yX)  ^max  (2“24) 

A  comparison  of  Equations  2-24  and  2-20  shows  that 
T  2 

-Ttiax  ^yy  Ttiax  ”  ^max^  (2"25) 

The  expression  on  the  left-hand  side  of  Equation  2-25  is  the 
autopower  spectrum  of  the  estimated  reference  channel,  which  is  automati¬ 
cally  normalized  to  the  numerical  value  of  the  group  coherence  (|  X  I2) 

'  max'  /* 
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A  comparison  of  Equations  2-23  and  2-25  shows  that 


n 

=xy 


h 

—max 


X 


max 


(2-26) 


The  expression  on  the  left-hand  side  of  Equation  2-26  is  the 
crosspower  spectrum  between  the  reference  channel  and  the  estimated 
reference  channel. 


E.  SUMMARY  OF  MINIMUM-POWER  ARRAY  PROCESSING 


The  salient  features  of  this  section  are  summarized  below: 

•  Equations  2-9  and  2-10  can  be  iteratively 

solved,  using  the  power  method,  to  give 

the  group  coherence  IX  I  2  and  the 

'  max 1 

MCrW;  the  MC!Xax  Sen'rate» 
the  maximum -coherent  reference 

channel  from  group  x 

•  MCF  li  ,  which  estimates  the  reference 

lllaX 

channel  by  linearly  operating  on  group  £, 
can  be  obtained  *rom  £Tnax  by  Equation  2-7 

•  The  minimum-power  array  processor  is 
the  MCF 

h 

—max 

^max 

and  the  power  in  the  corresponding  error 
channel  is  given  by 


min ' 


-  1 


-  1^ 


max ' 
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F.  THE  COMPLETE  SET  OF  m  EIGENVALUES 


Thus  far,  only  the  group  coherence  which  is  numerically 
equal  to  the  largest  eigenvalue  of  Equations  2-9  and  2-10  has  been  dis¬ 
cussed.  However,  the  generalized  eigenvalue  equations  have  m-1  other 
solutions. 

The  properties  of  the  entire  set  of  m  eigenvalues  and  the 
associated  eigenvectors  are  investigated  in  Appendix  A. 
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SECTION  III 
DATA  PREPARATION 


The  minimum-power  array  processing  technique  consists  of 
the  following  three  stages: 

•  Computation  of  the  cross  power  matrix  at  a 
chosen  frequency 

•  Computation  of  the  maximum  group- coherence 
MCF's  for  the  two  partitioned  sets 

•  Computation  of  wavenumber  responses  of  the 
two  group-coherence  MCF's  and  of  the  mini¬ 
mum-power  array  processor 

Presented  in  this  report  are  the  results  from  two  sets  of 
models  and  the  TFO  long-noise  sample.  1  The  TFO  crossarray  (Figure  III- 1 ) 
was  partitioned  into  four  sets  (Figure  III- 2)  and  group  coherence  was  studied 
at  various  frequencies. 

The  three  models  designed  at  0.  22584  Hz,  0.  52696  Hz,  and 
0.  82808  Hz  in  model  set  1  are  close  approximations  to  the  TFO  noise  field. 
Figure  III-3  shows  these  wavenumber  spectrum  models,  using  the  following 
conventions. 

A  solid  disk  represents  isotropic  energy  propagating  above  a 
certain  minimum  speed,  a  cylinder  represents  isotropic  energy  propagating 
at  a  fixed  speed,  and  an  arrow  represents  directional  energy  such  as  that 
generated  by  storms  and  earthquakes.  Background  mantle  P-wave  energy, 
due  to  its  high  apparent  speed,  is  modeled  by  an  8-km/sec  disk.  The  isotro¬ 
pic  fundamental  Rayleigh-wave  energy  is  modeled  by  the  3-km/sec  cylinder 
and  higher-order  modes  by  the  3-km/sec  disk.  Highly  directional  P  waves 
and  Rayleigh  waves  are  represented  by  the  arrows. 
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Figure  III  - 1 .  Array  at  Tonto  Forest  Seismological  Observatory  (TFO) 
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Figure  III-2.  Partitioned  Sets  of  TFO  Array- 
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Figure  IH-3.  TFO  Noise  Model  Set  1  in  Wavenumber  Space 
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Model  set  2  is  the  same  as  Model  set  1  except  that  1 -percent 
white  noise  is  added. 

The  crosspower  matrices  for  these  models  are  derived  from 
formulas  published  in  the  final  report  on  Seismometer  Array  and  Data  Pro¬ 
cessing  Systems.4  These  formulas  are  reproduced  in  Appendix  B. 

The  crosspower  matrix  for  TFO  data  is  computed  from  spectral 
estimates  of  the  TFO  noise  sample5  obtained  by  Bartlett- smoothing  and 
Fourier-transforming  the  auto-  and  crosscorrelation  functions  discussed  in 
an  Array  Research  semiannual  technical  report.  ^ 
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SECTION  IV 


DISCUSSION  OF  RESULTS 

With  respect  to  the  original  purpose  of  determining  seismometer 
inequalization,  the  study  of  group  coherence  filters  has  been  unsuccessful  in  achiev¬ 
ing  useful  results.  However,  the  following  detailed  discussion  of  results  does  illu¬ 
minate  the  properties  of  group  coherence  filters  and  indicates  that  they  will  be  use¬ 
ful  for  uncovering  and  analyzing  coherent  energy  between  subarrays. 

A.  PARTITIONS  A  AND  B 

Table  IV-1  shows  group  coherences  obtained  from  the  four  partitioned 
sets  A,  B,  C,  and  D.  Partitions  A  and  B  give  high  group  coherences. 

Minimum -powe r  array  MCF's  computed  from  TFO  data  at  0.  52696  Hz, 
corresponding  to  partitions  A  and  B,  are  plotted  as  vectors  in  Figure  IV-1.  An 
arbitrary  scale  factor  and  an  arbitrary  phase  reference  are  used  for  plotting  these 
vectors.  The  scale  factors  are  different  for  the  two  partitions. 

Both  vector  diagrams  show  that  the  MCF's  with  the  largest  weights 
are  located  near  the  center  seismometer,  which  is  not  included  in  these  partitions. 
Thus,  both  groups  x  and  seem  to  be  trying  to  predict  the  output  of  the  missing 
center  seismometer. 

Figure  IV-2  shows  wavenumber  responses  of  the  minimum -power 

arrays  (shown  in  Figure  IV-1  as  vector  diagrams).  These  wavenumber  responses  are 

computed  along  the  two  TFO  arms,  i.e.  ,  S37°W-N37°E  and  S53°W-N53°E.  Also 

shown  in  Figure  IV-2  are  the  inverted  K-line  wavenumber  power  spectra,*  which  are 

projections  of  2 -dimensional  power -density  spectra  onto  the  two  arms  of  the  TFO 
1  _ 

crossarray;  they  serve  as  references  for  checking  the  performance  of  the  minimum - 
power  arrays. 

*K-line  wavenumber  spectra  are  obtained  from  the  maximum  entropy  spectral 
analysis,  a  technique  developed  by  John  P.  Burg  and  presented  in  November  1967 
at  the  37tk  annual  SEG  meeting  in  Oklahoma  City.  The  properties  of  K-line 
wavenumber  spectra  have  been  discussed  and  extensively  illuminated  in  Array 
Research  Special  Report  No.  23.  * 
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The  broad  low  region  in  the  wavenumber  power  responses  of 
the  minimum-power  arrays  suggests  that  a  reasonable  noise  wavenumber 
power  spectra  should  have  a  corresponding  broad  peak.  This  is  not,  in  fact, 
the  case,  as  can  be  seen  by  comparing  these  wavenumber  responses  with  the 
inverted  K-line  wavenumber  power  spectra.  However,  a  basic  similarity 
between  the  two  shapes  can  be  seen,  as  both  the  minimum-power  array 
wavenumber  response  and  the  inverted  K-line  wavenumber  spectra  indicate 
that  the  ambient  seismic-noise  energy  sharply  drops  outside  the  3-km/sec 
dashed  lines. 

B.  PARTITIONS  C  AND  D 

Partitions  C  and  D  have  been  selected  to  provide  higher  reso¬ 
lution  in  the  wavenumber  response  than  achieved  by  partitions  A  and  B.  The 
minimum  separation  distance  between  the  two  groups  of  seismometers  in 
partitions  C  and  D  is  nearly  4  km,  which  is  two  times  the  minimum  separa¬ 
tion  for  partitions  A  and  B.  With  increased  separation,  the  responses  of  the 
two  group- cohe rence  filters  can  be  similar  only  over  a  much  narrower  wave- 
number  region  tnan  in  the  cases  of  partitions  A  and  B;  also,  as  a  result  of 
increased  separation  between  the  groups,  the  coherences  for  partitions  C 
and  D  are  low  (Table  IV-1). 
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:or  Diagram  of  Minimum- Power  Array  Processor, 
>  Noise  Data:  Frequency,  0.  52696  Hz 
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Figures  IV-3  and  IV-4  show  the  wavenumber  power  responses 
along  the  two  TFO  arms.  Minimum-power  array  wavenumber  responses  in 
model  set  1  generally  resemble  the  energy  distribution  in  the  model  but  do 
not  show  the  very  strong  directional  Rayleigh  wave  coming  from  N50°E 
(0.22584  Hz).  Minimum-power  array  wavenumber  responses  for  model  set  2 
show  poor  resolution.  Minimum-power  array  wavenumber  responses  obtained 
from  TFO  data  generally  resemble  the  corresponding  inverted  K-line  wave- 
number  spectra  obtained  from  the  same  data,  but  the  former  show  poor  reso¬ 
lution  and  are  unreliable  for  estimating  a  reasonable  wavenumber  power 
spectra.  This  last  conclusion  is  based  on  the  presence  of  extraneous  low 
regions  in  the  wavenumber  responses  at  0.52696  Hz  (TFO  data,  partition  C). 
Since  a  similar  effect  is  observed  for  synthetic  data  (model  set  2,  0.  52696Hz, 
partition  D),  this  phenomenon  appears  to  be  a  property  of  the  technique  and 
not  of  inequalization. 

C.  2-DIMENSIONAL  WAVENUMBER  POWER  RESPONSE 

Figure  IV-5  shows  2-dimensional  wavenumber  responses  of 
the  minimum-power  array  computed  at  the  three  frequencies  for  partition  D, 
model  set  1 . 

Comparisons  of  these  wavenumber  responses  with  the  wave- 
number  power  distributions  in  the  actual  model  (Figure  III-  1)  again  show 
that  the  technique  lacks  the  resolution  and  reliability  needed  for  estimating 
with  any  reasonable  accuracy  the  wavenumber  power  spectra;  at  least,  this 
is  true  for  the  array  groups  which  have  little  separation,  such  as  TFO. 

D.  GROUP-COHERENCE  MULTICHANNEL  FILTERS 

The  wavenumber  responses  of  group-coherence  MCF's  show 
very  poor  resolution  and,  since  no  conclusions  can  be  drawn  from  them,  they 
are  not  presented  in  this  report. 
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Comparisons  of  Minimum-Power  Array  Wavenumber  Responses  and 
Inverted  K-Line  Wavenumber  Spectra,  S37°W -N37°E  Arm 
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APPENDIX  A 


and 


THE  m  SOLUTIONS  OF  THE  GENERALIZED 
EIGENVALUE  PROBLEM 

The  generalized  eigenvalue  problem  is  stated  as 


G  n  G  =  I 

=  =xx  =  = 


(A- 1) 


where 


gt  Tn  n  _1  n  1  g  =  |\ 

=  [=xy  ~yy  =yxJ  =  = 


(A-2) 


1^  denotes  m  x  m  identity  matrix 
2 

|\|  denotes  m  x  m  diagonal  matrix  of  eigenvalues 
G  denotes  m  x  m  nonsingular  matrix  of  eigenvectors 

Section  II  of  this  report  shows  that  a  pair  of  MCF's  g^^  and 

h  exist,  corresponding  to  the  eigenvalue  |\  |2,  which  produce  two 

max  max 

maximum  coherent  channels  when  linearly  combined  with  groups  x  and 
respectively.  Similarly,  a  pair  of  MCF's  g^  and  lr  exist,  corresponding  to 
each  of  the  m  eigenvalues  of  the  set 

{|\.  i2  i  =  1,  .  .  .  ,mj 

where  g^  is  called  the  reference  MCF  and  lr  is  called  the  prediction  MCF. 

The  MCF  h.  is  related  to  the  MCF  by  Equation  2-7,  which  is 


h.  =  n 

-i  =yv 


Q  g. 
— yx  -“l 


(A-3) 


A- 1 
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Since  the  setjj^  |  i  =  1 . m|  is  the  m  x  m  matrix  G, 

the  set  H  =  jh.  I  i  =  l,...,m|is  annxm  matrix  given  by 

=  l—i  * 


8  ■  fiyy  =yx  S 


(A-4) 


To  complete  the  relationships  between  the  matrices  G,  H,  and 

i  ,2  . 

|  X  |  ,  the  following  two  equations  are  given: 


ht  n  h  =  |  x  | 2 

=  =yy  = 


(A-5) 


GT  Qx^  H  =  |  X  | 2 


(A-6) 


These  are  equivalent  to  Equations  2-25  and  2-26  for  the  m 


eigenvalue  case. 


•  Lemma  1 


The  m  reference  MCF's,  which  are  the  column  vectors  of  the 


m  x  m  matrix  G,  are  linearly  independent. 


Since  G  is  nonsingular,  its  column  vectors  are  linearly  in¬ 


dependent. 


Q.E.D. 


•  Lemma  2 


The  m  prediction  MCF's  which  are  the  column  vectors  of  the 

2 

n  x  m  matrix  H  are  linearly  independent  if  the  matrix  ]  X|  is  nonsingular. 

Rank  is  defined  as  the  maximum  number  of  linearly  independent 
column  or  row  vectors  in  a  matrix.  Since  all  the  columns  or  rows  of  a  non¬ 
singular  matrix  are  linearly  independent,  its  rank  equals  the  number  of  rows 

i  .2 

or  columns.  Thus,  the  rank  of  the  m  x  m  matrix  ]X  |  is  m  and  that  of  the 
n  x  n  matrix  fi  is  n. 

=yy 


The  ranks  of  the  matrices  on  the  two  sides  of  Equation  A- 5 

are  given  by 

p(hT  Syy  h)  =  p  (m2)  .  m  (A- 7) 

where  p  denotes  the  rank.  According  to  a  theorem  in  matrix  theory,  the 
rank  of  a  product  of  matrices  cannot  exceed  the  rank  of  any  of  the  component 
matrices  and,  since  p  (H)  =  p  (HT),  Equation  A-7  takes  the  form 


gives 


min 


2  m 


(A-8) 


Since  H  is  an  n  x  m  matrix,  another  theorem  in  matrix  theory 


p  (H)  £  min  (m,  n) 


(A-9) 


It  is  seen  from  Inequalities  A-8  and  A-9  that 


p  (H)  =  m  S  n  (A- 10) 

Equation  A- 10  says  that  the  number  of  linearly  independent 
vectors  in  H  is  m  and  that  this  number  cannot  exceed  the  dimension  (n)  of 
the  space. 

Q.E.D. 


•  Lemma  3 

The  column  vectors  of  the  m  x  m  matrix  G,  when  used  as 
reference  MCF's  on  group  x,  give  m  linearly  independent  outputs. 


Letting  the  outputs  of  the  MCF's 

random  variables  x.  and  x., 

-i  ■ 


x. 


(A- 1 1) 


and 


(A-12) 
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To  prove  the  lemma,  the  covariance  a. .  between  x.  and  x. 
will  be  shown  to  be  0  for  all  i  a  1,  ....  m  and  j  *  i:  J 


ij 


T 

=  x.  x. 

“i  ~1 


-  5)  (ij  2) 


T  T 
%  *  * 


=  h  Qxx 


(A- 13) 


T 

Equation  A-l  shows  that  ^  n  £.  =  0  if  i  ±  j. 

Thus,  it  can  be  seen  that  =  0.  Since  this  is  true  for  all 
1  =  1,  .  . .  ,  m  and  j  4  i,  the  theorem  is  proved.  The  term  "linear  independence' 
is  used  in  Lemmas  3,  4,  and  5  in  a  statistical  sense. 


Q.E.D. 


•  Lemma  4 

The  column  vectors  of  the  n  x  m  matrix  H,  when  used  as  pre¬ 
diction  MCF's  on  group  give  m  linearly  independent  outputs. 

This  lemma  can  be  proved  by  following  the  same  procedure  as 

T*  o 

used  in  Lemma  3  and  recalling  the  H  Q  H  =  |  X  |Z  (Equation  A-5). 

Q.  E.D. 

•  Lemma  5 


An  error  channel  is  defined  as  a  random  variable  obtained 
from  the  difference  between  the  outputs  of  the  reference  MCF  and  the  pre 
diction  MCF.  The  error  channels  generated  by  m  such  MCF  pairs  are 
linearly  independent. 
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To  prove  the  lemma,  the  covariance  0..  between  the 

ij 

channels  and  6 .  will  be  shown  to  be  0  for  all  i  =  1,  ....  m  and  j  ^ 


CT.  . 

1J 


G  . 
1 


^  i 

J 

i _ 

i 

X 

Gil 

-h. 

o 

-h. 

—l 

.  -yx 

=yyj 

“J 

°ij  (^i  -XX  — i  =yx) 

-  (gT  Q  -  hT  n  )  h. 

=xy  -i  =yy/  -j 


From  Equations  A-l  and  A- 5, 


T 

h 


n 

=xx 


,  T 

h.  n  h.  =  o 
-1  =yy  -j 


H  be  two 


(A- 14) 


(A- 15) 

error 


(A- 16) 


(A- 17) 
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Therefore,  Equation  A- 17  simplifies  to 


a. . 

ij 


0 

=yx 


x  T 

+  *i 


n 

=xy 


Equation  A-6  and  its  transpose  conjugate  show  that 


T 


h. 

“J 


n 

=yx 


0 


Substituting  this  result  in  Equation  A- 18  yields  a  =  0 

ij 

true  for  all  i  =  1, .  . .  ,  m  and  j  £  i. 

Q.  E.  D. 


(A-  18) 


.  This 
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APPENDIX  F 

FORMULAS  FOR  CROSSPOWER  SI  ECTRA  USED  IN  MODELS 


•  Crosspower  spectrum  due  to  directional 
energy:  ,  , 

-j  2  TT  f  •■krT  COS  0 
*12(f)  =  e  >V' 


•  Crosspower  spectrum  due  to  isotropic 
energy  propagating  with  fixed  speed 
(cylinder  model): 


where  represents  the  zero-order  Bessel 
function. 


•  Crosspower  spectrum  due  to  isotropic 

energy  propagating  with  speeds  above  V  : 

© 


* 


12 


(f)  * 


Seismometer  2 


|x|  *  Separation  between 
seismometers 

|  V  |  ~  Propagation  speed 

9  =  Angle  between  line  joining 
seismometers  and 
propagation  direction 

f  =  frequency 

Figure  B-l 


where  represents  the  first-order  Bessel  function. 

NOTE 


The  crosspower  spectrum  in  the  first  case  depends  on  the  seis¬ 
mometer  separation  |x  |  as  well  as  on  the  seismometer  orienta- 
tation  0  relative  to  the  direction  of  propagation.  In  the  latter 
two  cases,  however,  the  crosspower  spectra  are  independent  of 
the  orientation  0  and  depend  only  on  the  separation  |x|. 
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SUMMARY  OF  NOTATION 


Conjugate  transpose 
Contained  in 
Absolute  value 

Complex  multivariate  random  column  vector  of  m  elements 
|x.  |  i  =  1,.  .  .  ,  m|  which  are  complex  Fourier  transforms  of  the 
outputs  of  m  seismometers  in  group  x  at  fixed  frequency  f 

Complex  multivariate  random  column  vector  of  n  elements 
{Vj  |  j  =  1  *  •  •  •  *  n|  corresponding  to  group  jjr 

Complex  column  vector  of  m  elements  {g.  |  i  =  1,  .  . .  ,mf  which  is  the 
reference  MCF  associated  with  group  x 

Complex  column  vector  of  n  elements  |h.  |  j  =  1, . . .  ,n}  which  is  the 
prediction  MCF  associated  with  g~oup  y 

Error  channel,  which  is  the  difference  between  the  outputs  of  MCF  £ 
and  that  of  MC  F  ji 

)  Mean-square-error 

m  x  m  autopower  matrix  for  group  x 

n  x  n  autopower  matrix  for  group  ^ 

m  x  n  crosspower  matrix  for  groups  x  and  £ 
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n  x  n  identity  matrix 


I 

I 


G 

H 


m  x  m  identity  matrix 

•  2 
m  x  m  diagonal  matrix  of  eigenvalues  |\  | 

m  x  m  matrix  formed  by  group-coherence  MCF'sj^  |  i  =  1.  .  . 
n  x  m  matrix  formed  by  group-coherence  MCF'sjh.  |  j  =  1,  .  .  . 


p(H)  rank  of  matrix  H 

min(a,  b)  smaller  of  the  two  quantities  a  and  b 


_ UNCLASSI 


Security  Clmsification 


DOCUMENT  CONTROL  DATA  .  R  &  D 

r«£r  a  ^ —  — »  sssl  «sl  as  — » a  a  ^ , 

Texas  Instruments  Incorporated 

Science  Services  Division  Unclassified  _  _ 

P.  O.  Box  5621,  Dallas,  Texas  75222  26' °ROUP _ 

*  REPORT  TITLE  - - — — — — _ 

MINIMUM-POWER  ARRAY  PROCESSING  OF  THE  TFO  LONG-NOISE  SAMpt  it 
ADVANCED  ARRAY  RESEARCH  -  SPECIAL  REPORT  NO  12  LE 


«  OESCRIPTIVE  NOTE*  (7yP.  ol  r.porf  mnd  lnclu.lv .  dmt», ,  - - - - - - - - 

Special 

•  AUTHOR'S)  (Pint  nmno.  mlddlt  Inlllml,  /«•(  name)  - - — - - ■—  _ _ _ 

Burg,  John  P. 

Alam,  Aftab 

1 5  February  1 968  7*’  T0TAL  no4£fpaoe* F  N°  «« 

*».  CONTRACT  OR  GRANT  NO  - - — - -  I _ 

ORIGINATOR1*  REPORT  NUMDER1S) 

Contract  F33657-67-C-0708-P001 

b.  PROJECT  NO.  _ _ 

VELA  T/7701 

,e-  no<*’  (Any  num6#" m-y  »•  *••<*••<* 


OIITNI BUTION  1  1 

II.  *UPPLEVEN T ARY  NOTE*  .  ,  ’ 

12.  SPONSORING  MILITARY  ACTIVITY 

ARPA  Order  No.  624  Advanced  Research  Projects  Agency 

ARPA  Program  Code  No.  7F10  o{  ®?fe.n.se 


urogram  Code  No.  7F10  TK»P.‘V - w 

„  A*.T3)!(ly _ _ [  The  Pentagon,  Washington,  D.  C.  20301 

Vl  TW.  Teport  investigates  the  effectiveness  of  the  minimum-power  arrav  proceed 
tion.  th  ”  inequalizations.  The  technique  involves  parti- 

.o  ?haf  array  ‘nt0  tw°  groups  and  designing  MCF-s  for  each  group 

at  the  mean- square -error  between  the  two  MCF  outputs  is  a  minimum  unH^r  th* 
constraint  that  the  output  power  of  one  ?f  the  MCF' s  is  unity  The  two  MCF 

tfuTminirnum  8roup“coherence  filters:  the  difference  between  these  sets  is  known  as 
the  minimum-power  array  processor. 

Estimates  of  the  noise  wavenumber  spectrum  from  the  wavenumber  responses  of 

Z  J  *f  ,rVenr  /ureSP°n8e“  0fthe  two  group-coherence  MCF-s  are  very  similar 
(  eg.,  at  the  peaks  of  the  noise  wavenumber  power  spectrum)  7 

Seismometer  inequalization  was  to  be  determined  from  the  adjustment  in  weivh, 

and  phase  required  for  each  filter  so  that  the  wavenumber  responses  of  L  „-„ 

coherence  MCF-s  would  agree  with  a  reasonable  no.se  wavenumber  spctuTTow 

ever,  results  from  die  TFO  long-noise  sample  and  two  synthetic  models  shTw  , ha,  tte 

ma»mum  coherent  channels. 

conclusion  is  at  least  true  fc  small  arravs  such  as^FCh  lne<lualizatiGn.  This  latter 

DD  »  nov  es  1473  - - - - — ' 


Security  Classification 


UNCLASSIFIED 

Security  Classification 


